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Abstract. We study the nature of the dynamics in a first-order quantum phase transition between 
spherical and prolate-deformed nuclear shapes. Classical and quantum analyses reveal a change in 
the system from a chaotic Henon-Heiles behavior on the spherical side into a pronounced regular 
dynamics on the deformed side. Both order and chaos persist in the coexistence region and their 
interplay reflects the Landau potential landscape and the impact of collective rotations. 
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Shape-Phase transitions in nuclei are an example of quantum phase transitions (QPTs) 
in a mesoscopic (finite) system. These are qualitative changes in the properties of the 
system induced by a variation of parameters A in the quantum Hamiltonian H{X) [1]. 
Such ground-state transformations have become a topic of great interest in different 
branches of physics [2]. The competing interactions that drive these transitions, can 
affect dramatically the nature of the dynamics and, in some cases, lead to an intricate 
interplay of order and chaos. In the present contribution we study this effect [3, 4] 
in relation to first-order QPTs between spherical and axially-deformed nuclei [5], as 
encountered in the Nd-Sm-Gd region. We employ the interacting boson model (IBM) [6] 
involving {s, d) bosons with angular momentum L = 0, 2, representing valence nucleon 
pairs. The model has been widely used in describing QPTs [7] and chaos [8] in nuclei. 

In studying QPTs, it is convenient to resolve the IBM Hamiltonian into two parts, 
H = Hijii + Hco\ [9]. The intrinsic part (Ant) determines the potential surface V(/3,7), 
while the collective part (i?coi) is composed of kinetic terms which do not affect the 
shape of y(/3,7). Here (/3,7) are quadrupole shape parameters whose values (/3eq,7eq) 
at the global minimum of V (/3 , 7) define the equilibrium shape for a given Hamiltonian. 
Focusing on first-order QPTs between stable spherical (/3eq = 0) and prolate-deformed 
(/3eq > 0, 7eq = 0) shapes, the intrinsic Hamiltonian reads 

HUp)/h2 = 2{l-p^l5i)h,{nd-l)+P^Rlip)-R2iP): (1) 
H{i{^)/h2 = ^pIPo + pI-Pi, (2) 

where h2 = h2/N{N— 1). Here is the J-boson number operator, R2^{p) = y/ls'^ dj^ + 
pVl{dU^)^,^\ Pi=d^ ■ Jf - /3o2(/)2, P^^ = V2[5osUl + ^{dU^)f. The parameters 

that control the QPT are p and t,, with < p <^q^ and > 0, while /3o is a constant. 
H[^^{P) and hI^^{^) are the intrinsic Hamiltonians in the spherical and deformed phases, 
respectively. They coincide at the critical point pc = [5q ^ and ^c = 0: Hf^^ (pc) = H(^^ {^c)- 



The classical limit is obtained through the use of coherent states, rescaling and taking 
N ^ oo^ with \/N playing the role of h [8]. The derived classical Hamiltonian involves 
complicated expressions of shape variables (j8,7), Euler angles and their conjugate 
momenta. Setting the latter to zero, yields the following classical (Landau) potentials 

J8o2j82-pj8o2y2^i83cos37+i(l-i82)i84 , (3) 

+i[2(l-/3o^)+^(l+/3o¥]r + ^i8o^ (4) 

The (/3, 7) variables can be parametrized by Cartesian coordinates x = (5 cos 7 and y = 
j8 sin 7. The potential V\p) [V^'i^)] has a global spherical [deformed] minimum with, 
respectively, j8eq = [j8eq = \/2j8o(l + jSq^)"^/^ > 0, 7eq = 0]. At the spinodal point (p*), 
V^(p) develops an additional local deformed minimum, and the two minima cross and 
become degenerate at the critical point pc (or ^c)- The spherical minimum turns local in 
V^'{^) for ^ > and disappears at the anti-spinodal point (^**). The order parameter 
j8e^ is a double- valued function in the coexistence region (in-between p* and £,**) and 
a step-function outside it. The potentials V(j8, 7= 0) = V{x,y = 0) for several values of 
^ , p, are shown at the bottom row of Fig. 1 . The height of the barrier at the critical point 
is Vfc = /z2 [1 — ( 1+ Pq) ^/^]^/2. Henceforth, we set j8o = 1.35, resulting in a high barrier 
Vi,/h2 — 0.231 (compared to Vi,/h2 — 0.00lS in previous works [8]). 

The classical dynamics of L=0 vibrations, governed by flint, can be depicted conve- 
niently via Poincare sections. These are shown in Fig. 1 for selected energies and control 
parameters. For p =0, the system is integrable, with V\p =0) o<: jSqjS^ + ^(1 — /3o)/3'^. 
The sections for p =0.03 in Fig. 1, show the phase space portrait typical of an an- 
harmonic (quartic) oscillator with two major regular islands, weakly perturbed by the 
small pcos37term. For small /3, V^{p) ^ j8^— p-\/2/3^cos37. The derived phase-space 
portrait, shown for p =0.2 in Fig. 1, is similar to the Henon-Heiles system (HH) [10] 
with regularity at low energy [panels (bi)-(b2)] and marked onset of chaos at higher 
energies [panels (b3)-(b5)]. The chaotic component of the dynamics increases with p 
and maximizes at the spinodal point p* = 0.546. The dynamics changes profoundly in 
the coexistence region, shown for p = 0.65, 0.741 and =0.05 in Fig. 1. As the lo- 
cal deformed minimum develops, robustly regular dynamics attached to it appears. The 
trajectories form a single island and remain regular at energies well above the barrier 
height Vh, clearly separated from the surrounding chaotic environment. As ^ increases, 
the spherical minimum becomes shallower, the HH-like dynamics diminishes and disap- 
pears at the anti-spinodal point ^** =0.354. Regular motion prevails for ^>^**, where 
the section landscape changes from a single to several regular islands and the dynamics 
is sensitive to local degeneracies of normal-modes [4]. 

The quantum manifestations of the rich classical dynamics can be studied via Peres 
lattices {xi.Ej} [11]. Here are the energies of eigenstates \i) of the Hamiltonian and 
Xi = ^y2{i\hd\i)/N. The lattices can distinguish regular from irregular states by means 
of ordered patterns and disordered meshes of points, respectively. The particular choice 
of Xi can associate the states with a given region in phase space through the classical- 
quantum correspondence j8 =jC'f^jc/ [3]. The Peres lattices for L=0 eigenstates of flint 
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FIGURE 1. Poincare sections (upper five rows) depicting the classical dynamics of Hj^^j(p) (1) and 
H'^^{^ ) (2) with 112 = 1, Po— 1-35, for several values of p < p,. and ^ > ^c- The bottom row displays 
the corresponding classical potentials V (p ) (3) and V ) (4) and the five energies, below the domain 
boundary ^^^^^(jS = \/2, 7), at which the sections were calculated. The Peres lattices {xi,Ei}, portraying 
the quantum dynamics for eigenstates with L = and = 80, are overlayed on the classical potentials 
V{x,y = 0). They exhibit sequences of regular states in the vicinity of the deformed well, consisting of 
bandhead states of the ground ^(/T = 0), I3"{K = 0), l3"y^{K = 0), I3'Y{K = 0) bands, etc. 



with = 80, are shown on the bottom row of Fig. 1 , overlayed on the classical potentials 
V{x,y = 0). For p =0, the Hamiltonian (1) has U(5) dynamical symmetry with a solvable 
spectrum Ei = 2^2 [PqN— 1 + ( 1 — /3q For large A^, the corresponding Peres lattice 

coincides with V'{p =0),a trend seen in Fig. 1 . Whenever a deformed minimum occurs 
in the potential, the Peres lattices exhibit regular sequences of states, localized in the 
region of the deformed well and persisting to energies well above the barrier. They 
are related to the regular islands in the Poincare sections and are well separated from 
the remaining states, which form disordered (chaotic) meshes of points at high energy. 
The regular L = states form bandheads of ^ = rotational bands. Additional ^-bands 
corresponding to multiple /3 and 7 vibrations, can also be identified. An example of such 
regular K = 0,2 bands for H^^i at the critical point, is shown in Fig. 2(a). The states in 
each band have nearly equal values of (n^/), indicating a common intrinsic structure. 
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FIGURE 2. Peres lattices {xi,Ei} for A^ = 50,L = 0,2,3,4 eigenstates of //^^(p = p,.) =H"t{£, = ^c), 
Eqs. (l)-(2), with /12 = 1, j3o = 1-35 [panel (a)] and additional collective terms involving 0(3), 0(5) and 
0(6) rotations [panels (b), (c), and (d)]. The classical potential shown, is the same in all cases. Notice in 
panels (a)-(b)-(c), the well-developed rotational bands (K — Q, L — 0,2,4) and {K = 2, L = 2,3,4) formed 
by the regular states in the deformed phase, which are distorted in panel (d). 

The collective part of the Hamiltonian (Hcoi), which does not affect y(/3,7), is 
composed of the two-body parts of the Casimir operators of the groups 0(3), 0(5) 
and 0(6) [9]. These orthogonal rotations are associated with the Euler angles, 7 and 
/3 degrees of freedom, respectively. Fig. 2 shows the Peres lattices corresponding to 
L = 0, 2, 3, 4 eigenstates of H^jii at the critical-point, plus added rotational terms one at a 
time. As seen in Figs. 2(b)-2(c), the 0(3) and 0(5) terms preserve the ordered ^-bands 
of H[jii, Fig. 2(a). In contrast, the regular band-structure is strongly disrupted by the 0(6) 
term [Fig. 2(d)]. The latter couples the deformed and spherical configurations [12] and 
mixes strongly the regular and irregular states. These results demonstrate the advantage 
of using the resolution H = Hint + ^coh since a strong 0(6) term in the collective part 
can obscure the simple patterns of the dynamics disclosed by the intrinsic part. 
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